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Abstract. The Casimir energy of a dilute dielectric cylinder, with the same light- velocity 

) (where 

si + e 2 J 

Q ■ £i, £2 are the dielectric constants of the cylinder and of its environment), and numerically to 
higher orders in £ 2 . The first part is carried out using addition formulas for Bessel functions, 
■ and no Debye expansions are required. 
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^ '. 1 Introduction 

Zero-point fluctuations of quantum fields give rise to forces, which are regarded as manifes- 
tations of the Casimir effect (for reviews, see e.g. refs.0). From the theoretical viewpoint, 
one of the most daunting aspects of the evaluation of Casimir energies, even for highly sym- 
metrical boundaries, is its sheer difficulty. Many mathematical methods has been developed, 
but even the simplest ones demand considerable efforts. 

At the core of several of these techniques one finds uniform asymptotic expansions 
— also called Debye expansions — of Bessel or Ricatti-Bessel functions appearing in integrals 
over momentum-like variables. This fruitful method was used as early as — at least — the 
time of ref. 0], and has been repeatedly revisited in a huge number of articles, often in 
the framework of other regularization schemes (see e.g. ref. || and refs. therein). However 
reliable, the whole Debye expansion technique is a time-consuming process and the search for 
computational alternatives might be of interest |4]] . This is, precisely, one of the motivations 
of the present letter. Our purpose is to take further the exploitation of summation theorems 
for Bessel functions, started in ref. Q for cases with spherical surfaces, and apply it to a 
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problem with a cylindrical boundary. 

We are considering a material cylinder of radius a, infinitely long, placed along the z- 
axis, with permitivitty and permeability e±, fi±, surrounded by a medium with permitivitty 
and permeability 62,^2- For such surfaces, a special situation is the case where the light 
velocities in both media — interior (1) and exterior (2) — are the same, i.e., 

eiHi = £2^2 = c" 2 , (1) 

where c is the common light-velocity. Since any variation in e affects fi, this is called 
dielectric- diamagnetic opposed to the purely dielectric one, in which [i\ — ^2 — 1 

but the velocity has to change. Dielectric-diamagnetic conditions are often desirable as they 
cause the frequency equations to simplify and some divergences to cancel out. In a QCD 
context, e and /x refer to colour permitivitty and permeability (see []5J and refs. therein). 
Illustrations of the dependence of the interquark potential on the boundary conditions for a 
string model have been provided in ref.|J. 

In refs. [[7), |§ and || the regularized Casimir energy per lateral unit-length for an 
infinite dielectric-diamagnetic cylinder has been studied. Up to the order of 

e = {^)\ (2) 

V£i + e 2 J 

the energy has been shown to vanish within all the tested degrees of numerical accuracy. 
The next contribution, which is of the order of £ 4 , has been found — to our knowledge, for 
the first time — in ref.||. 

In our minds, medium 2 will be pure vacuum and medium 1 a very tenuous dielectric, 
which means 62 = = 1 and S\ — 1 <C 1. As a result, the £ 2 parameter, defined by 
eq.(fj), is a small number. According to ref.0, the eigenfrequencies uj coming from the 
Maxwell equations for this problem are given by the zeros of some equations of the type 
f n (k z ,LU,a) = 0, n G Z — eqs.(2.3)-(2.5) in ref.[0] — . Further, in cases where the relation ([[]) 
holds, f n takes the form 



f n (k z ,u,a) = -a c X 



2 -2\e( £ i + £ 2) 2 



4:6x62 



eV 2 n (Xa) + 



V n (x) EE (J n H n )'( X ), (3) 



where £ 2 is given by (0), J n , H n are Bessel and Hankel functions, and 



uj = cJX 2 + kl (4) 
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Every A belongs to the eigenfrequency set of the projected two-dimensional problem — say 
A — , while — oo < k z < oo, i.e., the values of k z are continuous without any restriction. 
Before regularizing, the Casimir energy per unit-length (£c) is given by the mode sum 

£ C=^J2 I ^n,m,fc z , ^n,m,k z = C ^X^ + k 2 . (5) 

Z n,m J Z7r 

The n-index is the angular momentum number, while m describes the remaining degree of 
freedom, i.e., labels the different A- values at a given n. 

The present work is organized as follows. In sec. 2 we follow ref. and evaluate the 
energy density (in momentum space) up to the order of £ 2 , by a modified Bessel function 
summation theorem, and resorting to the properties of Meijer G functions. Then, we show 
that the integration of yields a vanishing result. Sec. 3 is devoted to an alternative 
approach based on a zeta function prescription for the initial mode sum like in refs.0, |i"0 |. 
Apart from proving to be easier, this technique paves the way to the numerical calculation 
of higher order contributions. Our conclusions are given in sec. 4. 



2 Density method 

We begin by reviewing the procedure used in ref. @ and obtaining an expression for the 
Casimir energy. The mode sum is first represented, as usual in these cases, by a contour 
integral 
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fn,as(k z , LO) 



(6) 



where the integration contour C consists of a straight line parallel to, and just to the right 
of, the imaginary axis, (—zoo, +zoo) closed by a semicircle of an infinitely large radius in the 
right half-plane. The branch line of the function y?(A) = 



A 2 + k\ is chosen to run between 
i\k z \ and i\k z \ on the imaginary axis. In terms of y = ImA we have 



y 2 _ 


U2 


A 2 


-y 2 , 



i\/y 2 - k 2 z , 



y > h, 
\y\ < k z , 
y < -k z . 



(7) 



Noting that the argument of the logarithm is an even function of iy, 

he 



reduces to 



E, 
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E 



27r n"oo J 
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(8) 
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where we have expressed (|3|) explicitly on the imaginary axis. Integrating with respect to k z 
we obtain the Casimir energy per unit length (Ec) as an integral over y, namely 

£ c = ^ E / dyy 2 d y ln[l-e(ydy(I n (ay)K n (ay))) 2 ] . (9) 

47F n=-oo J ° 

In the following subsections, we use the approach of ref. || to evaluate this expression. 
2.1 Density calculation to the order of £ 2 

Having integrated k z out, the integrand in expression (Q) may be interpreted as the density 
of Casimir energy with respect to the parameter y = iX. This density may be evaluated by 
expanding in terms of £ 2 , the first contribution being 

g V\y) = ^- dyy 2 d y [yd y (I n (ay)K n (ay))} 2 . (10) 



71= — CO 



We shall now show that g^ 2 ' can be calculated explicitly by a variant of the method shown 
in ref. [|J. Specifying the identity 8.530.2 of ref. [11] to Hankel solutions Z n = H^> = H n , 
and choosing the v parameter equal to zero, we obtain the summation theorem 



H (mR(p,r,(p)) = J2 J n( m p) H n {mr) e mip , R(p, r, <p) = p 2 + r 2 - 2pr cos <p. (11) 

n=— 00 

Performing the change m — > im, and selecting the special case p = r, it becomes 



K (mR(r,ip)) = £ I n (mr) K n (mr) e inip , R(r, <p) = ry/2(l - cos <p) = 2r |sin (ip/2)\ . 

n=— 00 

(12) 

Differentiating with respect to m, using the property K' (z) = —Ki(z) together with the fact 

that K n = K_ n , and setting m = 1 afterwards, we have 

00 

-R(r,<p)K 1 (R{r,<p))= ]T r (/„ K n )\r) e^. (13) 

n=— 00 

Recalling the orthogonality of the imaginary exponential functions, we arrive at 

— / dp[R(r,< P )K 1 (R(r,< P ))] 2 = £ [r(J n i^ n )'(r)] 2 = F(r), (14) 

In order to proceed, we rename r into ay, and do the variable change u = |sin (</?/2)|. After 
differentiation with respect to y, we realize that the sum in (|i~0|) is given by 

0( 2 ) ( y ) = ?£L J y 2 [^(/ n (a?/)K n (a?/))] 2 = ^y- / du^=^=y 2 d y (ayK 1 (2ayu)) 2 , 

(15) 



Thus, we have turned the problem of calculating an infinite angular-momentum summation 
into the evaluation of a definite integral of a transcendental function. Next, by writing the 
product of Bessel functions appearing in Eq.(^) in terms of the Meijer G function ||12| , gW 
can be evaluated explicitly. (See Appendix for the definition and some simple properties of 
this function.) First, we write 



u'aVK^ayu) = uW ^2 Q; -1, 0, 1; Aa 2 y 2 u 2 ) = ^G™ (J; 0, 1, 2; Aa 2 y 2 u 2 ) 

(16) 

where, in the last step, we have made use of the identity fl48|). Differentiating ( |T6"D with 
respect to y and substituting in (|15|) we obtain 

<J t2) (v) = -^I'^tS ?" (i;0,o,i;4a 2 yV) . (17) 

Using eq.5.5.2(5) in [O, and some straightforward manipulation one gets: 



2.2 Calculation of the energy to order £ 2 . 

We now turn to the question of deriving the Casimir energy. There are two possibilities: 

1. To perform the u integration in flTBD after integrating over the y variable. Unfortunately, 
this turns out to be divergent. Considering the integral 

— J y 2 dy dud y [y sin uK^y sin u)] (19) 

and interchanging the order of integration, one arrives at 

roc ^ 1 

/ dyy 2 d y [ysinuK l {2ysinu)] = . (20) 

Jo 12 sin u 

If we now try to do the u-integration, the integral diverges. This shows that some further 
regularization is called for. In fact, in sec. ^ we will go through the same sort of calculation, 
but with the advantage of having applied zeta function regularization from the outset. 

2. Direct integration of g <y2 \y)- We use the following identity from ref. |T2[ (Vol 1, page 
215). 

POO 

/ dyy- a K v {2jy) G™ (a h . . . , a p ; b u . . . , b q ; xy) 

J „ 



\ r m,n+2 ( _ V_ \ ( 21 ) 

2 < - T p+2,g \ a 2 ' a 2 ' ' ' ' ' ° P ' X ' ' ' ' ' *?' X 



[2 _ 

In order to take advantage of this formula, we note that Ki(y) = W — e y . Changing to a 

2 V 

2 2 

variable t = (Jr ^ and inserting K i , we can cast the energy per unit-length into the form 

(22) 

Although we have inserted Ki , as a mere technicality to help calculate the energy, one can 
think of using it as an exponential regulator^ (see ref.[|J). However the convergence of the 
integral shows that the density we have derived is already regularized in some sense. One 
can now use ( pip to get: 

^ (2) = ~J~2~~/ = ^44 (0, 1 1, 2; 3 - 3 - , §, 1; 4) . (23) 
4a z Tr-s/ircr r V 2 2 2 2 2 cr r . / 

In order to check the asymptotics as o r — > 0, we use the property (fE|), together with the 
asymptotics (again from ref.fO) G(x) = Odxl' 3 ) as x — > 0, for p < q, and /3 = maxjRe^} 
for h = 1, . . . ,m. In our case we simply have 

£ c (2) oc lim — 0{\a 2 T \) = lim 0(a r ) = 0. (24) 



Thus, the £ 2 -term is shown to vanish, confirming the conclusions of refs. 0, || and 0, 
without recourse to numerical evaluations. 



3 Complete zeta function regularization 

In this section we will take a different approach, based on the application of the complete 
zeta function method (see e.g. refs.0, |Hj) to the initial mode sum (M). The use of zeta 



functions for regularizing such sort of sums dates from the time of refs. [13]. In the version 



we shall now apply, the regularized value of the Casimir energy per unit-length is 

£c = lim -hc(n(D=3)(s) (25) 
s— ►— l 2 

where the zeta function (n(D=3) for the whole set of tu-modes in the three-dimensional problem 
— say fl — is given by 

^)W=Er|(^f)"' (26) 



3 Actually, it is not difficult to show that applying an exponential regulator in the form e arUJ to 
before carrying the k z integration, yields the same result as the one we derive. 
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First, one assumes that s is large enough for this function to make sense, with the final aim 
of setting s = — 1 at the end (usually, one introduces in (|25| ) an arbitrary mass scale, but, in 
this problem, it turns out to be unnecessary). Taking into account (|5|), we write 



/oo rlh r 
IT- K,m + k 
, -oo ^vr 



Let's consider the zeta function for the projected two-dimensional problem, i.e., for the A 
eigenmode set: 



CA(D=2) i 



(J 



CO 

E 

n=0 



do = 1, 

d n = 2, for n > 1, 



(2f 



where Cn(c) stands for the nth partial-wave zeta function 

oo 



(29) 



m=l 



Bearing this in mind, we put (^) as 
Cn(D= 3 )(s) = — B ^) Ca(d= 2 )(s - 1) 



CA( 7+i ^ + ( ln ( 2 ) " 5) Ca(d= 2) (-2) + Ca(W-2) + + 1) 



(30) 



1 

27 



where an expansion around s = — 1 has taken place. This was the method applied in ref . [|T0| . 

Eqs. (|28|), (|29|) hold for Re cr > 1, but they will have to be analytically continued to 
the neighbourhood of a = —2 (note that a = s — 1). Such an analytic continuation is carried 
out by the contour integration method of refs.|T3|, |3[. To begin with, one takes 



a 



a a ( n (v) = 77— / duu a 1 ln [fn(uj] , for Re cr > 1 , 



(31) 



where f n (u) = f n (k z ,u,a) with u = Xa, and C is a circuit in the complex u-plane enclosing 
all the positive zeros of f n (u). In the desired limit this contour will be semicircular, with 
the straight parts along the imaginary axis, and adequately avoiding the origin. The first 
step (see e.g.flnj) is to examine the asymptotic behaviour / n>as (w) of f n (u) for \u\ — > 00. If 
fn,as( u ) has no roots inside of C, we leave the eq.(|3T|) unchanged by setting 

fn(u) 



a- u U(j) = — I duu-*- 1 ln 
2m Jc 



(32) 
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Going back to ([|), we observe that, for large x, V\{x) — O (x 4 ). Then, one can write 



fn(u) = /n.asO) 

Therefore, eq.(^) translates into 

a^CJa) = — [ duu^ 1 In 
2ni Jc 



fn >aS (u) 



{£l + E2? C " 2 u\ (33) 



El E2 VT 2 a 4 



(34) 



After realizing that only the vertical parts of C — where u = e ±tn ^ 2 y — are actually con- 
tributing to the integration, eq.(|34"D yields 



Cn((T)=a ^ sin ^j ^ dyy-*- Hn{l-e V (4 K n )'{y) 



for -1 < Re a < 0, 



(35) 



where we have used V^(±iy) 



(I n K n ) (y) , being I n , K n the corresponding modified 
Bessel functions (note that this y is dimensionless) . All this has validity near a = — 1, but 
we still need some further work in order to reach the neighbourhood of a = — 2. 



3.1 Calculation to the order of £ 2 

Let £c = E £c P ^ £ 2? \ an d analogously for the involved zeta functions. Then, 

p>i 



poo 

Ag>\*) = / dyy— 1 [y (I n K n )'(y)f p , for -1< Re a < 0. 
Jo 



(36) 



Note that we are commuting a ^-expansion with a process of analytic extension which 
sidesteps cr-poles (i.e., s-poles). Yet, since the ^-dependence has no problematic traits, 
this should be correct, and we write 



CA(D=2)(a) = £dnC» =ECi(l 2 )(^K 2P , 



n=0 



p>l 



A2 P ) 



1 a <r . fna\ 
--a — sin — 
p it \ 2 J 



J2 d nA^(a). 



(37) 



n=0 



If we just want to keep the terms ~ £ 2 in Sc, it will be enough to maintain the p = 1 
contribution, which can be rewritten in the way 



.(2) 

>A(D=2) ^ 



C.a//» o.i'Tj = -<i' ~ Sill 



77 V 2~ 



/ dyy^~ 1 F(y), F(y)= £ [y(W(y)] 2 , (38) 
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where we have also taken into account Q2"8D and the fact that £_ n (er) = Cni®)- An integral 
representation of the F(y) is already available in eq. (fH]) . From there, we proceed as in the 
derivation of eq. (|15|) , i.e., we do the variable change u = | sin(<£>/2)| and find 

j 2 f 1 du 



. i 8y f 1 du 9 ~, 

Hy) = ^r - 7 ==u 2 K 2 1 (2uy). 

7T JO \/l — U 



Vl-u 

With this, we go back to eq. and focus on the integral 

rl du 



dyy-°- 1 F(y) = - 



,r I dyy-° +1 Kl{2uy). 
o 



(39) 



(40) 



7T JO y/\ — U 2 

The ^-integration is evaluated with the help of formula 6.576.4 in ref.[|lTJ. Then, the remain 
ing w-integral is immediate using forumla 3.251.1 in the same book. As a result, 



r 




r 2 




r( 


:-s) 


r(2-a)r 




) 



(41) 



which has a zero of order one at a = —2 by virtue of the singularity of T ( £ ^) • Putting it 
into eq. (38) and expanding near a = —2, we find 



a + 2) 2 + 0((a + 2) : 



(42) 



which provides the desired analytic extension to Re a = —2. The crucial point is that it has 



(2) 

a zero of order two at a = — 2 and, therefore, Ca(d=2)( — 2) = and £ 
together with eqs. fl25|) and (0), leads to S c = + O (£ 4 ), i.e., 



a(d=2) (-2) = . This, 



(2) 

c 



(43) 



which was numerically found in ref.0 (see also [§, |]). 



3.2 Higher-order corrections in £ 2 

In order to know new corrections in £ 2 , one has to keep the next £ 2 -terms in eqs. (0), (j37|). 
However, unlike Affl(a), the A^ p \a) integrals with p > 2 are already finite at a = —2, 
because [y(I n K n )'(y)] 2p ~ ^ as y — > oo. Thus, the restriction to Re a > —1 in (|36|) is 

4 Incidentally, the comments made by one the authors of the present letter led to the correct numerical 
figure in ref. jf) . 
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caused only by the presence of p = 1, while, for p > 2, it suffices to numerically evaluate all 
the necessary A^(<t)'8 at a = -2, i.e., 

1 oo 

-2) = --^£^4 2p) (-2). 



A2p) ' 
t 'A(D=2) 



(44) 



P a n=0 

A posteriori, we have verified that each term decreases quickly enough with n and that the 
n-summation has a numerically acceptable behaviour. Then, the finiteness of each of these 

'—2) = 0, and the ^^-contributions to Sc are 



u(2p) 



sums confirms that Ca(_d=2)^ 



1 

2 
-Tic 



.(2p) 

n(D= 
An pa? 



1 



(-1)^ 

oo 



Rc-clf2. 2) '(-2)«* 



for p >2, 



(45) 



n=0 



where the meaning of Cnm=3)( s ) * s obvious. When p = 2, including all the n- values up to 



120, we have found ^ d n A<£\-2) ~ 0.19108. This and formula © yield 

he 



n>0 



4 4) e 4 



-0.0076028 ^£ 4 , 



(46) 



in agreement with ref. 0. As remarked there, the negative sign means that the involved 
Casimir forces are attractive. Physical implications concerning the flux tube model for 
confinement have been discussed in that work. 

In fact, the higher p, the fewer terms are needed in the n-series for obtaining reliable 
figures. For p > 2 we have found many contributions, but we just list the first ones: 



p 


3 


4 


5 


6 


7 




£g p) a 2 /(hc) 


-0.0022637 


-0.0010807 


-0.0006202 


-0.0003972 


-0.0002737 





As argued in refs.0 or ||, the special value £ 2 = 1 should reproduce the perfectly- 
conducting case Sc ( px .) a 2 /(he) = —0.01356 . . . 0. Taking all the contributions up to 
p = 7, we obtain Ec a 2 / (he) ~ —0.01224, with a 10% relative error. This is not too surprising, 
as the ^-expansion comes from a logarithimc series, and a slow numerical convergence at 
£ 2 = 1 is expectable. Including all the terms up to p = 200 we have found Sc a 2 / (he) ~ 
—0.01354, with a 0.15% relative error. 



4 Conclusions 

The ultimate consequences of any result about Casimir effect are not easy to foresee, as the 
domain of applicability of this concept has been expanding beyond what could be considered 



10 



'traditional' areas of field theory. For instance, we have recent examples of these ideas in 
spacetime evolution and quantum cosmology |T7|]. Proposals haven even been made about 



possible ways of extracting work from the vacuum energy ||18|| . 

In the present letter we have confirmed the expectation that the ^-contribution to 
the Casimir energy for a dilute-dielectric cylinder, infinitely long, and under the condition 
of light-velocity conservation, would have to vanish. Numerically speaking, this had been 
noticed with very high accuracy in several articles, starting with ref. j7|, but in the present 
letter we have been able to derive it as an exact result (eqs. fl24|) and fl43|) ). Another new 
aspect lies in applying the method developed in [||, i.e., the use of summation theorems for 
infinite series of Bessel functions, which has spared us the handling of Debye expansions (see 
also the application of this method in ref . J05] and , but this time in connection with the 
problem of refs.[pl|l). 



Moreover, by a numerical evaluation, and within the complete zeta function regu- 
larization framework, we have reobtained the ^-contribution calculated in ref. |J, (our eq. 
([46])), which is negative. This constitutes the first deviation from zero and shows that, 
although at a higher order, the Casimir energy of this system would tend to contract the 
cylinder. Even higher corrections in £ 2 have also been found (table in sec. [3]). 

Our spectral zeta function has been constructed like in ref.[|l(|. Other variants of 
the zeta function procedure, which differ from ours at some particular steps, are also in 
circulation (e.g. ref.JTH] or p|). We regard them as slightly different formulations of one 



common underlying principle. In particular, ref.||15|| illustrates the advantages of dealing 
with the total zeta function as a whole object, rather than a series of partial-wave zeta 
functions. 

A Appendix: The Meijer G function 

Here we state some facts about the Meijer G function, which is defined by the integral 

a-(a, u x)- 1 Ids nr =1 rfe- 3 )n^r(i- aj + 3 ) 

g„ (*., (.,», x)-—j L ds n^ra-tj + .jnjl^rfo-.) 1 • (47) 



The different integration paths L can be found, for example, in [12], as most of the other 



properties we use. By simple variable changes one may prove numerous identities such as: 

x n G™(a liat7 6urt, x) = G™(a Usi + n, b hst + n, x) (48 
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and 

G™(a m , b nst , ~) = G£ n (l - 6 list , 1 - a hst , x) (49) 



which we have used in sees. 2.1,2.2. 
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